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Abstract 

So far fluid mechanical Nambu brackets have mainly been given on an 
intuitive basis. Alternatively an algorithmic construction of such a bracket 
for the two-dimensional vorticity equation is presented here. Starting 
from the Lie-Poisson form and its algebraic properties it is shown how 
the Nambu representation can be explicitly constructed as the continuum 
limit from the structure preserving Zeitlin discretization. 



1 Introduction 

The Hamiltonian description has long been a matter of interest in fluid mechan- 
ics and it has enabled theoretical insight to turbulence [1, 2] and also practical 
applications [3, 4]. In a series of articles a generalization of Hamiltonian fluid 
mechanics in terms of Nambu mechanics [5] has been proposed [6, 7, 8, 9, 10]. 
Perspectives for the numerical analysis of fluid motion based on this formalism 
have also been given [7, 11, 12]. Not much is known however about the math- 
ematical background of the Nambu brackets in fluid mechanics and of generic 
PDEs, specifically a general algorithm for finding these brackets is lacking as 
mentioned in [7] . Addressing this topic we present here a constructive derivation 
of the Nambu bracket for the example of the two-dimensional barotropic vortic- 
ity equation. Moreover we discuss the generalized Jacobi identity, the validity of 
which has not been investigated thoroughly so far for fluid-mechanical Nambu 
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brackets, or has been completely ignored. The overall result may hence be of 
interest given the increasing number of publications on Nambu representation 
in fluid mechanics. 

The further part of section 1 gives an overview of the state of the art, com- 
prising the Zeitlin truncation scheme for the vorticity equation [13], Nambu 
brackets [5], and a general method by Bialynicki-Birula and Morrison [14] for 
their construction in the discrete case. Using the latter in section 2 a Nambu 
representation of the vorticity equation is constructively derived from the un- 
derlying algebraic properties. Based on our findings we conclude with discussing 
the generalized Jacobi identity. 



1.1 The Zeitlin truncation scheme for the vorticity equa- 
tion 

Wc consider the motion of an ideal, incompressible fluid on a compact domain 
CI. In Lagrangian view, the phase space is the cotangential bundle of the Lie 
group of volume preserving smooth diffeomorphisms on fl. This description can 
be reduced to an Eulerian one with phase space given by the dual of the Lie 
algebra sdifF(f2) [15]. In two dimensions the resulting dynamical system may be 
expressed through the barotropic vorticity equation 

% = -JiM. (1) 

72 



Here ^ is the vorticity of the fluid and tp the stream function related by C = V"'V' 

da db da dh 

dxi dx2 dx2 dxi 



The Jacobi operator is defined by J(a, b) := -S^-^ — and the velocity of 

^ .7^7/ rjx^ OX'} OXr> OX\ 



the fluid is given by v = '5x1^' ^'^^ simplicity all functions are assumed 

to be smooth. 

Under suitable (e. g. periodic) boundary conditions (1) can be written in 
Lie-Poisson form 

§ = {C,H}po.. (2) 
Here the Lie-Poisson bracket of the barotropic vorticity equation is defined by 

{^i,^2}poE:=y^^dACj(^^,^) (3) 

for differentiable functionals J-\ , Ti of the vorticity. The Hamiltonian is given 
by the kinetic energy of the fluid 

n:=\\ dAv2 = 1 f dA(VV^)2 = -1 /■ dAC^. (4) 

In the following we consider a rectangular domain of length 27r in each direction 
with periodic boundary conditions, i.e. the torus T^. Using the fact that the 
Lie algebra sdiff(T^) can be approximated by su(n) [16, 17, 18], a spectral 
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truncation scheme for (1) preserving the Lie~Poisson structure was proposed by 
Zeithn [13]: 



dCi _ n ^ 1 . (2-K 
'dt ~ ^2^ 



^ lc2 ^ ^) (i e /„). (5) 

kG/„ ^ ^ 

Here Ck := /t^ dA C(x)e~''' '' denote the Fourier coefficients of the vortic- 
ity. The wave vectors accounted for are located on the integer grid without the 
origin 

|i = (^l,^2) eZ2\{0} - ^- < < ^- j , (6) 

bearing a rectangular cutoff with — 1 modes in total (n odd) . The norm and 
the vector product of the indices arc determined by :— kf + fc| and i x k := 
11^2 — 12^1 respectively. The vertical dash in (5) stands for the dimension-wise 
modulo function onto the grid. 

The truncated form of the energy as a function of the Fourier coefficients is 

^^=^(27r)^E i^Ci^C-k. (7) 

ke/„ 

This approximates the continuum form of the energy (4) which follows using 
Parseval's identity Jt^ dA/(x)g(x) = X^kez^ /kffk and the reality condi- 
tion Ck = C-k- It can easily be checked that the truncation (5) is of Lie-Poisson 
form ^ 

^ = {Ci,i^}oo.. (8) 
The bracket of this discrete dynamical system is given by 



iJ,kG/„ ^Ci 9Cj 



^ . ,'27:. .W dF,dF2 



for scalar- valued diffcrentiable functions Fi, F2 of the vorticity modes. This Lie- 
Poisson bracket is completely determined by the structure constants (Q!y)ij.ke/„ 
of the underlying Lie algebra su{n) : 

where S denotes the Kroncckcr delta. In the continuum limit n ^ 00 they 



converge to the structure constants 



rk 
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r(ixj)<5i+j,k (11) 



'J ■ (27r)2 

of sdiff(T2). yielding the continuum vorticity equation (1) in spectral form. 
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1.2 Nambu brackets and their algebraic construction 

Recall that the Nambu form of a general dynamical system with prognostic 
variable z reads [5] 

= {2, 7^1,^2}. (12) 

In consistence with [7] we define the Nambu bracket {•,-,•} by the properties 
trilinearity, total antisymmetry, and the Leibniz rule; Hi and T-ii denote gener- 
alized Hamiltonians. Moreover a Nambu bracket defines an algebra of conserved 
quantities if and only if it fulfills the generalized Jacobi identity [19], see sub- 
section 2.3. 

In [6], a Nambu bracket for the vorticity equation has been given: 

(.r,..r.,^3}™:^/„..f.(f.f)^ (13) 

Inserting the enstrophy Casimir 

£:=^JjAC' (14) 

as third argument recovers the Lie-Poisson bracket (3) (assuming suitable bound- 
ary conditions): 

{J-l, J-2, S}pBE ~ {J'l, J'2}pDE- (15) 

Yet, in [6] neither was a claim made about the generalized Jacobi identity nor 
was an algorithm for deducing the Nambu bracket (13) given. 

For ODEs, such an explicit construction method from the starting point of a 
Lie-Poisson bracket was presented in Bialynicki-Birula and Morrison [14]: Let 
{cx^j) be the structure constants of a semi-simple n-dimensional Lie algebra g. 
If Zk {k = 1, . . . ,n) denote coordinates on the dual space g*, the Lie-Poisson 
bracket is given by 

{FuF,}= a%zu^^^ {F,,F,eC-'{g*)). (16) 

The Killing form Kij := '^ik'^ji is invertible with inverse denoted by K^L 

A Nambu bracket can then be defined through 

{F,,F2,F3}:= N,,k^^^ (Fi,F2,F3GC°°(0*)) (17) 

. ^ ^ dzi dzj dzk 

with Nambu tensor 

n 

iV,,fe :=^a^^.if,fc. (18) 

1=1 
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By definition a Casimir C of a Poisson bracket satisfies {C, •} = 0. The classical 
example in the semi-simple Lie-Poisson case is the quadratic Casimir 



1 

C := - ^ K^^z.z,. (19) 

Inserting it as third argument into the Nambu bracket reproduces the Lie- 
Poisson bracket: 

{Fi,F2,C} = {Fi,F2}. (20) 

Note that the total antisymmetry of the Nambu tensor and therefore also of the 
bracket follows from the antisymmetry and the Jacobi identity of the structure 
constants. 



2 Construction of a Nambu bracket for the dis- 
crete and continuum vorticity equation 

Unfortunately the method discussed above cannot be used directly for construct- 
ing Nambu brackets for PDEs, since the Killing form is generally not defined 
for infinite-dimensional Lie algebras. Specifically the Killing form 

k,lGZ2\{0} ^ ' keZ2\{0} 

corresponding to the structure constants (11) of the spectral vorticity equation 
diverges. The difficulties in the application of this method to hydro dynamical 
brackets have also been noted by [7]. Therefore in the following we make use 
of the structure preserving truncation by Zeitlin in order to derive the Nambu 
bracket of the vorticity equation according to the method of Bialynicki-Birula 
and Morrison. This comprises the computation of the Casimir and the bracket 
as well as the verification of their convergence in the continuum limit. 

2.1 The quadratic Casimir for the Zeithn discretization 

We compute the Killing form corresponding to the structure constants (10) of 
su(n): 



k,ieA, 

1 71^ ^ v-^ , / 47r \ , , s , X 

= -2(^'^(i+J)l".o 2^(1-^°%— '^""j) (ije/„). (22) 



ke/„ 

By the orthogonality relation 



cos f—k ■ l \ = nHi,o - 1 (23) 
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and since the origin is not part of the grid, we get X^ke/ (iT^ k) = — 1. 
Thus the Kilhng form and its inverse read 

1 n** .. (27r)'^ 

= -2 (2^'^('+j)l"'° ^'^"^ = -2^^%+j)|„.o- (24) 

Therefore, up to the normahzing constant r the Casimir (19) is 

equal to the truncated form E of the enstrophy 



= ^ E ^^^"CiCi = ^(27r)2 E CiC-i =^ (25) 
This indeed approximates the continuum enstrophy 

f dAC^ = i(2^)2 ^ CkC-k, (26) 

which shows that the algebraic Casimir property carries over to the conservation 
property of the dynamical PDE system. 

2.2 The Nambu bracket for the Zeithn discretization and 
its continuum hmit 

According to (18) we define the scaled Nambu tensor 

^'jk - I E "y^'ik = - sin (^^i X <5(i+j+k)l„,o (27) 

with associated bracket 

p PI 1 ^ . f27r. \ dFidF2dF3 

for functions Fi,F2,F3 of the vorticity modes. By construction the antisym- 
metry properties follow from the antisymmetry and the Jacobi identity of the 
underlying Lie-Poisson bracket (9), which is recovered by inserting the trun- 
cated enstrophy: 

{i^l,i^2,i^}0DB = {^^1,^^2}0DE. (29) 

In combination with the truncated energy the corresponding dynamical system 

^ = {Ci,i/,£;}ooE (30) 

thus also represents the Zeitlin discretization, but in a form where the conser- 
vation property for the enstrophy is obvious. Of course there are n — 2 more 
independent Casimirs in the Zeitlin discretization, but the goal here is to identify 
the relation to the continuum Nambu bracket (13) which employs the enstrophy. 
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The above discussion carries over straight to the continuum case: Assume 
the functions Fi on the n-modal system converge to the functionals J^i on the 
oo-modal system, which by definition are equal to the functionals J-i of the 
vorticity field: 

lim F,((Ck)ke/J = ^«((Ck)kGZ^\{o}) = -^.(C) = 1, 2, 3). (31) 

n— >oo 

Then the here defined Nambu bracket (28) converges to the bracket (13) of the 
vorticity equation^ : 

lim{^^i,F2,F3}ooE = 2^ (1 X j)Ji+j+k,o ^ ^ ^ 

= {^1, ^2, -T^sjpDE- (32) 

This shows, that the Nambu vorticity bracket discussed in [6] can be constructed 
from the starting point of the Zeitlin discretization in Lie-Poisson form along 
the ideas stated in [14]. Moreover, this discussion reveals the algebraic origin of 
its characteristic properties. This includes the antisymmetry property and the 
fact that it reduces to the Lie-Poisson bracket when used with the enstrophy. 

2.3 The generalized Jacobi identity 

A generalization of the Jacobi identity has been included in the definition of 
Nambu brackets by [19]. While the classical form of this identity ensures that 
the Poisson bracket of two conserved quantities is again a conserved quantity, the 
generalization does the same for a Nambu bracket of three conserved quantities. 
For fluid mechanical systems however this issue has not yet been thoroughly 
discussed. Here as well as in (application-oriented) publications, e. g. [5, 6, 7, 
10, 11], the generalized Jacobi identity is not included in the definition of a 
Nambu bracket. Now we assess its validity in the case of the vorticity equation. 
The generalized Jacobi identity for a Nambu bracket {■, •, •} reads 

{{i^i, F2, Fa}, F4, F5} + {F3, {Fi, F2, F4}, F5} (33) 
+{F3, F4, {Fi, F2, F5}} ~ {Fi, F2, {F3, F4, F5}} = 0, 

which must hold for arbitrary functionals Fi, F2, F3, F4, F5. For a totally anti- 
symmetric Nambu tensor N defined through 

rp p f;-V/V ^P^dF^dF, 

{F4,F2,F3} _ ^7V.,, — — — , (34) 
the generalized Jacobi identity (33) takes the form 

W.'=f^ifit)9F4aF5 aF3W,^itl|lt)9F5 




dzi dzp dzq dzi dzp dz, 



g 



^It is used that the derivatives are related through (^^T~—u: = ,-"*".,> (* = 1,2,3). 
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dzi dzp dzq dzi dzj dzk I 

Terms involving second-order derivatives of F^^ F^^ F^ cancel and, if N is inde- 
pendent of z, by the total antisymmetry of N this equation reduces to 

d^Fi dF2 dFi d^F2 \ dFs dF^ dF. 



dzidzi dzj dzi dzidzj J dzk dzp dzq 



^ = 0. (36) 



Thus in the case of a constant Nambu tensor, the generalized Jacobi identity 
(33) is equivalent to the identity 

N.,jkNlpq + N.jqNlkp + N,,pNlqk = (37) 

for all i,j,k,l,p,q in the finite or infinite integer range. 

However, as remarked in [20], the generalized Jacobi identity is not neces- 
sarily fulfilled for the Nambu tensor constructed with the method of Bialynicki- 
Birula and Morrison. A simple counterexample for the discrete bracket (28) as 
well as for the continuum bracket (32) is given here: With the wave vectors 
i= (1,0), j = (0,1), k= (-1,-1), 1 = (1,0), p= (-1,1), q= (0,-1), we get 
in both cases for the first term in (37) a non-vanishing value and for the others 
zero. Consequently, none of these Nambu brackets fulfills the generalized Jacobi 
identity. 



3 Conclusion 

The presented article contains a discussion on the interpretation of the vorticity 
equation Nambu bracket in terms of its Lie-Poisson bracket. The algebraic 
Casimir of the Zeitlin discretization was computed explicitly and shown to be 
proportional to the truncated enstrophy. It was further demonstrated how the 
Nambu bracket with its characteristic properties arises naturally from the Lie- 
Poisson form of the vorticity equation. The existence of a Nambu representation 
for the vorticity equation could be traced back to the fact that a structure 
preserving truncation with semi-simple Lie algebra and enstrophy as classical 
Casimir exists. Without contradicting this result the generalized Jacobi identity 
was shown not to apply. 

It would be interesting to generalize the construction method for the use 
with vorticity Casimirs of higher degree instead of the enstrophy, with poten- 
tial perspectives in numerical applications. Furthermore, while for geophysical 
fluid dynamics the here treated two-dimensional vorticity equation is already of 
interest, a future discussion could also include more sophisticated fluid mechan- 
ical systems, such as e. g. the three-dimensional incompressible or shallow-water 
equations. However, since no structure preserving discretizations exist in these 
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cases [13], this is presumably more challenging. It may still be hoped that the 
presented method can shed light on an interesting facet of the formulation of 
fluid mechanics. 
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